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Abstract. We prove that minimal graphs (other than planes) are parabolic in 
the sense that any bounded harmonic function is determined by its boundary 
values. The proof relies on using the coupling introduced in [2] to show that 
Brownian motion on such a minimal graph almost surely strikes the boundary 
in finite time. 



1. Introduction 

By a minimal graph, we mean a complete minimal surface-with-boundary, the 
interior of which can be written as a graph over some open subset of the plane. 
(We assume that our minimal graphs are connected.) A famous result of Bernstein 
(see [3]) states that the only minimal graphs over the entire plane are planes (that 
is, the graphs of afhne functions). In this paper, we prove that all other minimal 
graphs (aside from the planes just mentioned) are parabolic, in the sense that any 
bounded harmonic function is determined by its boundary values. In particular, 
this means that the boundary of such a surface is non-empty. We obtain this 
result as a corollary of proving that Brownian motion on any non-planar minimal 
graph almost surely has finite lifetime, which means that it almost surely strikes the 
boundary in finite time. Our proof of this is essentially an application the coupling 
introduced in [2]. 

With this result, we settle a conjecture (in the affirmative) due to Meeks that 
any minimal graph over a proper subdomain of the plane is parabolic. Related 
questions have been studied by Lopez and Perez [1], who prove that a non-flat 
properly immersed minimal surface-with-boundary that lies above a sublinear graph 
and that has Gauss map contained in an open hyperbolic subset of the sphere is 
parabolic. 

2. A DESCRIPTION OF THE COUPLING 

As mentioned, our proof leans heavily on the results of [2j- We summarize the 
relevant results, specialized to the case of a single minimal surface-with-boundary, 
here; the interested reader is encouraged to consult the original paper for more 
details. 

Let M be a minimal surface-with-boundary. Then, for any point (xo,yo) in the 
interior of M x M, a coupled Brownian motion is a continuous stochastic process 



Date: October 3, 2008. 

2000 Mathematics Subject Classification. Primary 53A10; Secondary 58J65 60H30. 
Key words and phrases, minimal surface, minimal graph, parabolicity, Brownian motion, 
coupling. 

The author gratefully acknowledges support from an NSF Postdoctoral Research Fellowship. 

1 



2 



ROBERT W. NEEL 



(xt,yt), starting at (xo,yo), such that the xt and yt marginals are both Brownian 
motions on M, defined until the first time xt — yt or either of the marginals strikes 
the boundary of M. In order to describe the particular coupling we will use, we 
need to introduce some processes associated to any coupled Brownian motion. First, 
let r t = distR3 (x t , yt ) be the (extrinsic) distance between the particles. Next, let 
to : M — > § 2 be the Gauss map of M. The coupling will be governed by the relative 
positions of three unit vectors (which we think of as points in § 2 ), namely m(xt), 
m(y t ), and a{xt,yt) — (x< — yt)/\xt — yt\- We will refer to the position of these 
three unit vectors as the configuration of the system. 

In what follows, we will use W z to denote a Brownian motion (on the real line) 
with time parameter z. The particular Brownian motion may vary from use to 
use, and we will use W z to denote a second such Brownian motion (which need 
not be independent of the first). In [2], it was proved that there exists a coupled 
Brownian motion on M, starting from any (xo,yo)j such that the semi-martingale 
decomposition of r t can be written as 

dn = y^fdWt + ^-dt 

where / and g are non-negative functions of m{xt), m(yt), and a(xt,yt) (with 
one caveat, described below) satisfying several additional properties which we now 
describe. (That / and g depend only on the configuration was not required of the 
coupled Brownian motions considered in [2], but it's easy to see that this can be 
arranged.) First, the inequality / > g holds everywhere, and we take a moment to 
explain the significance of this. If we introduce the time-change r(t) = J Q fds, then 
the martingale part of r T is a Brownian motion (we somewhat abuse notation and 
use r T to denote what should be rj( r )) and the inequality / > g means that r T is 
dominated by a two-dimensional Bcssel process. Further, if the inequality is strict 
at any instant, then the domination is also instantaneously strict. 

Second, we can characterize the configurations where / = g. They all occur when 
m(xt), m(jjt), and a(x y ,yt) lie on the same great circle. Let S be the subset of 
(S 2 ) 3 such that all three points lie on a great circle in § 2 . On any simply-connected 
subset of S, we can define a coordinate for the great circle by the signed distance 
from a(x t ,yt) (obviously, there are two possible choices of coordinate related by a 
change of sign). Then let 9 be the coordinate at m(x t ) and let tp be the coordinate 
at m(yt). We would like 9 and ip to be continuous, and thus we will think of them 
as being defined only up to multiples of 2ir. Since we will only be concerned with 
the values of trig functions at 9 and tp, defining them only up to 2tt won't cause 
any problems. We have that, on this component of S, 

f = (sin# — Asinp) 2 and g = (cos 9 — Acosip) 2 
where A — sgn (cos(6* + ip)) . 

Here the choice of whether A is —1 or 1 at a point (x,y) with cos(6> + <p) = 
depends on the geometry of M x M near (x, y), and this is the one caveat to our 
statement that the evolution of r t (instantaneously) depends only on r t and the 
relative positions of m{xt), m(y t ), and a(xt,yt)- Nonetheless, the determination of 
A at such points isn't relevant to this paper. Continuing, we see that / = g — 
when 9 — ip and A = 1. Also, when cos(0 + ip) = 0, we see that there are 
two possibilities for / and <?, corresponding to A = 1 and A = — 1, and for each 
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possibility we have / = g > 0, with one exception. Namely, when cos(# + tp) = 
and 9 = ip, we have that A = 1 gives / = g = and A — ~1 gives f = g = 2. 

The above gives us good control of the evolution of r t for configurations in S. 
Our next goal is to extend this to a larger set of configurations. We introduce a 
sequence of subsets of (§ 2 ) 3 which we will use throughout the remainder of the 
paper. Let Si be the set where m(xt) and m(yt) are at least some fixed (small) 
positive distance ci apart in the § 2 metric. (We won't specify a value for ci; rather, 
we will assume that it satisfies various properties as we go.) Let £4 be the subset of 
S where cos(# + </?) = and where m(xt) and m(y t ) are at least distance 2ci apart. 
Let S3 be a closed neighborhood of S4. Let S2 be an closed set containing 53 on 
which I cos(# + ip)\ < C2 for some positive constant C2, and such that the boundary 
of S2 is positive distance from the boundary of S3, in the product metric on (§ 2 ) 3 . 
Further, we assume that S 2 and S3 are chosen so that 



and so that the boundaries of all four sets are a positive distance from each other. 
Much as we do for ci, we will feel free to adjust 53 and S2 as we go. 

We extend 9 and ip from £4 to be smooth functions on all of S2. (Note that this 
extension is not the same as used in [2], but this is a more convenient definition 
of 9 and (p for our present purpose.) On S3 we have that f — g > C3 cos(6> + <p)\ 
for some positive constant C3, while on Si \ S3, we have that / — g > C4 for some 
positive constant C4. 

Now that we have good control of / — g in terms of 9 + ip, we need to understand 
how the configuration evolves. Again, we will simply cite facts from [2]. Both 
m(xt) and m(yt) are time-changed spherical Brownian motions, with the time- 
changes given by the integrals of —K along the respective paths on M. Further, 
a(xt,yt) is a semi-martingale on S 2 , such that both the quadratic variation and 
the drift grow as a rate bounded from above by a multiple of 1/r 2 . This gives 
us control over how quickly the configuration can move from one region of (§ 2 ) 3 
to another in terms of K and r t . In the other direction, we will need to know 
that the configuration doesn't spend too much time in S3. Recall the time-change 
r(t) = J* * / dt. Since / is bounded above and below by positive constants on Si, so 
is dr I dt. In particular, if we assume that the configuration of our coupled Brownian 
motion never leaves Si , then a process has infinite lifetime in the original time t if 
and only if it also has infinite lifetime in the r time-scale. Also, all of our earlier 
estimates on the evolution of processes in the t time-scale hold, up to constants, 
in the r time-scale, assuming the configuration stays in Si. We have the following 
semi- martingale decomposition for (9 + p) T , valid on S2, 



where A takes only the values —1 and 1 and extends our earlier definition of A on S 
to S 2 . In particular, A is —1 on some open subset of S 2 , 1 on another open subset, 
and its value on their common boundary depends on the geometry of M x M. Here 
ei and €2 are bounded functions on S2, with common bound e > that can be 
made arbitrarily small by shrinking S 2 . Further, a and b are functions depending 
on the Gauss map at x and y with the property that both a 2 and b are bounded 
pointwise by a multiple of the sum of the absolute values of the Gauss curvatures at 
x and y, with the consequence that / a 2 dr and J bdr along any path (in M x M) 
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are bounded by a multiple of the sum of the integrals of the absolute values of the 
Gauss curvatures along the x t and y t marginals. 

3. The coupling on minimal graphs 

We now have enough information about the behavior of the coupled Brownian 
motion to prove our main theorem. Note that by embedded we mean injectively 
immersed. 

Theorem 1. Let M be an embedded, complete minimal surface-with-boundary 
which is not a plane. Further assume that, outside of a compact subset of M , 
the Gauss map of M is contained in an open, hyperbolic subset of the sphere. Then 
Brownian motion on M almost surely hits the boundary in finite time, and thus M 
is parabolic. 

Proof: The theorem is clear when M is (isometric to) a subset of the plane. 
Thus we assume that M is not flat, and we have that the Gauss curvature, K < 0, 
has only isolated zeroes. 

Let B t be a Brownian motion on M , started at a point p in the interior. Let a be 
the first hitting time of the boundary; we will always stop the Brownian motion at 
the boundary. We wish to prove that a is almost surely finite. The induced process 
on the Gauss sphere, m(B tAtT ), is a time-changed Brownian motion on § 2 , with the 
time-change given by u(t A cr) = / * A<T —KoB v dv, the integral of the absolute value 
of the Gauss curvature along the path. Our assumptions on the image of m imply 
that u(t A cr) almost surely converges to a finite limit as t — > oo and thus that the 
process m{B tAa ) almost surely converges. Further, using that K has only isolated 
zeroes, we see that the limiting distribution of the process on the sphere, which we 
denote m{B ooAa ) (despite the fact that Boo act is not itself well-defined a priori), 
does not charge any points. 

We now suppose that a is not almost surely finite and derive a contradiction. 
Under this assumption, we can find starting points for Brownian motion on M 
for which the probability of {cr = oo} is arbitrarily close to 1. Further, because 
fn{Boo/\a) does not charge points, among such starting points we can find two 
points xo and yo with the following further property: there exist subsets of S 2 , 
Ti and T 2 , such that the four sets ±ri and ±r 2 are positive distance apart and 
such that Brownian motion started at xo has normal vector that stays in Ti for its 
entire lifetime with probability arbitrarily close to 1, and similarly for yo and r 2 . 
It follows that we can find points xo and yo, with corresponding sets Ti and T 2 , 
such that the coupled Brownian motion described above, started from (xo, 2/0)1 has 
probability at least 1/2 of never hitting the boundary and of (m(x tA(T ), m(?/ tA(T )) 
staying in Ti x T 2 for all time. Further, because M is embedded, note that x t 
and y t cannot couple (that is, r t cannot strike zero) when {m{xt ArT ),m{yt Aa )) is in 

ri x r 2 . 

Next, we wish to study our adequate coupling, as described in the previous 
section, started at (xo,yo)- Let cr be the first time the process either hits the 
boundary or has (m(xt) , m(yt)) exit Ti x T 2 . Let £1 be the set of paths (under 
the adequate coupling started at (xo,yo)) for which a = 00. Then paths in ft have 
infinite lifetime, since they also never couple, and the probability of f2 is at least 
1/2. We can choose a set Si, as in the previous section, such that all paths in ft 
have configurations which stay in Si for all time. Since / is bounded from below 
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on Si and the paths in Q have infinite lifetime, it follows that r T almost surely 
has infinite quadratic variation on Q. Because r T is dominated (from above) by a 
two-dimensional Bessel process and is bounded from below by a Brownian motion, 
we know that infinite quadratic variation means that r T almost surely hits every 
positive level infinitely often. 

The contradiction we're working toward is to show that r T almost surely does 
not hit every positive level infinitely often on Q. To see this, we introduce the 
change of coordinate p = logr and the time-change s — f (1/r; 2 .) dr. Then we have 

(1) dp. = dW a - 1 ~ 2 g// da. 

Further, we have subsets S 2 , S3, and S4 of Si as described in the previous section 
such that 

d{9 + if) s = (2 + ei) dW s + ardW s + [A(2 + e 2 ) + br 2 ] ds 

on S 2 . Here we note that a simple computation, using that the integrals of the 
Gauss curvature along the two marginals are almost surely finite, shows that the 
integrals J a 2 r 2 ds and J br 2 ds are almost surely finite. 

Note that, because / is bounded from below on Si, our inequalities for g/ f imply 
similar inequalities for 1 — g / f . In particular, we see from Equation |T]) that p s 
is a Brownian motion with non-positive drift and that this drift is bounded from 
above by a negative constant on Si \ S3, We wish to show that this drift, in some 
average sense, spends significant time away from zero. We partition s-time into the 
intervals /„ = [n — l,n) for all non-negative integers n and consider the behavior 
of (9 + ip) s over these intervals, since this will allow us to estimate the drift of p s . 

We proceed by a series of special cases. First, assume that a and b are identically 
zero, and consider the process until a. We have d(8 + tp) s — (2 + ei) dW s + A{2 + 
e 2 ) ds. We wish to prove that the drift (of p s ) has a fixed probability of being less 
than a negative constant over I n . This follows from showing that, during /„, the 
configuration has a positive probability of spending some positive amount of time 
outside of S3. By a change of measure using Girsanov's theorem, we can eliminate 
the drift of (9 + ip) s so that it becomes a time-changed Brownian motion, with time 
change (relative to s) bounded below by (2 — e) 2 and above by (2 + e) 2 . Recall 
that I cos (8 + ^)| < C2 on S2. Thus we see that the configuration has a positive 
probability of leaving S 2 within time 81 > 0, where Si can be made as small as 
we wish by choosing C2 to be small. Next, note that if the configuration is in 
Si \ S2, it has a positive probability of taking at least time 5 2 > to return to S3. 
This follows from the fact that the boundaries of S 2 and S3 are positive distance 
apart along with our earlier observations about the evolution of the configuration, 
adjusted for the change to s-time. Finally, we need some understanding of how these 
estimates relate to one another. Our coupled Brownian motion and its associated 
configuration process are not necessarily Markov. Nonetheless, these estimates (on 
how quickly the configuration exits S 2 and how quickly it returns to S3) hold for any 
process, as long as the coefficients of its semi-martingale decomposition obey the 
right estimates, and are thus independent of the past. (Intuitively, these estimates 
are the result of implicitly comparing our processes with the "worst-case" processes, 
and these comparison processes are Markov.) We conclude that, after choosing our 
various constants and sets appropriately (by which we mean the Ci, the 5i, and the 
Si), there is a positive probability, call it 7, that the drift of p s will be less that —6, 
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for some positive S, on any interval /„. Further, by the same reasoning as above, 
the estimate holds on each interval independently of the others. Thus, there is a 
countable sequence X n of independent Bernoulli random variables, each of which 
is equal to 1 with probability 7 and with probability 1 — 7, such that 

(2) / -L-ilLd»<-8 £ X n 

where the integral is along paths and \z\ is the largest integer less than or equal 
to z. 

We now wish to lift the requirement that a and b are identically zero. Instead, 
we assume that Jj a 2 r 2 ds and Jj br 2 ds are less than some small, positive e over 
any interval I n with n < a. If e is small enough, then with probability arbitrarily 
close to one, the contribution of the terms involving a and b to the evolution of 
the configuration can be made arbitrarily small. Thus by assuming e is sufficiently 
small, and perhaps making 7 and S slightly smaller and further adjusting our con- 
stants Si and Ci and our sets Si, the estimate in Equation j2|) still holds. Applying 
the strong law of large numbers to the sequence X±,X2, ■ ■ ., we see that, almost 
surely for paths in ft, J —(1 — g/f)/2ds goes to negative infinity linearly. That is, 
there exists some positive constant c such that 

^ P 1 - ill A J_ 

(6) / ds + cs — > — 00 as Sq — > 00 



for almost every path in Q. 

Finally, we come to the general case, where we make no additional assumptions 
on a and b. As mentioned above, we know that the integrals J a 2 r 2 ds and J br 2 ds 
are almost surely finite. Thus, for any choice of e, we have that almost every path in 

has Jj a 2 r 2 ds and Jj br 2 ds bounded by e for all but finitely many n. Throwing 
away finitely many intervals doesn't effect the asymptotic behavior of the drift, so 
we see that Equation J3|) holds in general. 

The law of the iterated logarithm implies that | W S(j | almost surely grows sub- 
linearly, in the sense that W So /sq — > almost surely as sq — > 00. It follows from 
Equation ([I]) and Equation ([3]) that p s — * —00 as s — ► 00 almost surely on Q. This 
means that for these paths p s has a last time above 0, which in turn means that 
r T has a last time above 1. This contradicts our earlier observation that almost 
every path in Q has infinitely many excursions above any level, and we conclude 
that Brownian motion on M almost surely strikes the boundary in finite time. 

The last claim of the theorem, that this implies parabolicity, follows from the 
standard representation of bounded harmonic functions in terms of Brownian mo- 
tion evaluated at stopping times. □ 

Since any graph is embedded and has its Gauss map restricted to a hemisphere, 
any minimal graph, other than a plane, satisfies the assumptions of the above 
theorem. Thus we've succeeded in proving the following result. 

Corrolary 2. Any minimal graph which is not a plane is parabolic. 
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